Determining initial-state fluctuations from flow measurements in heavy-ion collisions 



(N 

o 

(N 
C 

in 



o 

=5 



> 
o 

o 



Rajeev S. Bhalerao, 1 Matthew Luzum, 2 and Jean- Yves Ollitrault 3 

1 Department of Theoretical Physics, TIFR, Horni Bhabha Road, Colaba, Murnbai 400 005, India 
2 CEA, IPhT, (CNRS, URA2306), Institut de physique theorique de Saclay, F-91191 Gif-sur-Yvette, France 
3 CNRS, URA2306, IPhT, Institut de physique theorique de Saclay, F-91191 Gif-sur-Yvette, France 

(Dated: January 20, 2013) 

We present a number of independent flow observables that can be measured using multiparticle 
azimuthal correlations in heavy-ion collisions. Some of these observables are already well known, 
such as t>2{2} and i>2{4}, but most are new — in particular, joint correlations between vi, V2 and 
«3. Taken together, these measurements will allow for a more precise determination of the medium 
properties than is currently possible. In particular, by taking ratios of these observables, we construct 
quantities which are less sensitive to the hydrodynamic response of the medium, and thus more 
directly characterize the initial-state fluctuations of the event shape, which may constrain models 
for early-time, non-equilibrium QCD dynamics. We present predictions for these ratios using two 
Monte-Carlo models, and compare to available data. 

PACS numbers: 25.75.Ld, 24.10.Nz 



INTRODUCTION 

Thermalization of the matter produced in ultrarela- 
tivistic nucleus-nucleus collisions results in strong col- 
lective motion. The clearest experimental signature of 
collective motion is obtained from azimuthal correlations 
between outgoing particles. It has been recently real- 
ized [l| that fluctuations due to the internal structure 
of colliding nuclei (previously studied in the context of 
elliptic flow |2j), followed by collective flow, naturally 
generate specific patterns which are observed in these az- 
imuthal correlations. In this Letter, we propose a number 
of independent flow measurements and study the possi- 
bility to constrain models of initial-state fluctuations di- 
rectly from these experimental data. 



FLOW OBSERVABLES 

Correlations between particles emitted in relativistic 
heavy-ion collisions at large relative pseudorapidity A?/ 
are now understood as coming from collective flow 
According to this picture, particles in a given event are 
emitted independently according to some azimuthal dis- 
tribution. The most general distribution can be written 
as a sum of Fourier components, 
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where v n is the n flow harmonic 4] and ^ n the corre- 
sponding reference angle, all of which fluctuate event-by- 
event. 

The largest flow harmonic is elliptic flow, V2 I5| .which 
has been extensively studied at SPS j™, RHIC [7H9( , and 
LHC [lfj. Next is triangular flow, 1)3 TQ, which together 
with t>2 is responsible for the ridge and shoulder struc- 



tures observed in two-particle correlations ll|,|l2j. In ad- 



dition, quadrangular flow,_U4, has been measured in cor- 
relation with elliptic flow [ljj E I- Finally, directed flow, 
can be uniquely separated 15( into a rapidity-odd 
part, which is the traditional directed flow [Taj, and 

a rapidity-even part created by initial fluctuations fl7j j . 
which has only been measured indirectly [l8j ]. 

In practice, one cannot exactly reconstruct the under- 
lying probability distribution from the finite sample of 
particles emitted in a given event. All known information 
about v n is inferred from azimuthal correlations. Gener- 
ally, a /c-particle correlation is of the type 

v{ni, ri2, ■ . ■ , n k } = (cos (m<pi + ■■■ + n k ip k )) , (2) 

where m, . . . ,n k are integers, tpi, . . . , ip k are azimuthal 
angles of particles belonging to the same event, and an- 
gular brackets denote average over multiplets of particles 
and events in a centrality class. Since the impact pa- 
rameter orientation is uncontrolled, the only measurable 
correlations have azimuthal symmetry: n\ + . . . + n k = 0. 

In this work, we are interested in the global event 
shape. The average in Eq. @ is thus taken over all mul- 
tiplets of particles. More differential analyses (i.e., re- 
stricting one or several particles to a specific p t interval) 
are left for future work. The average in Eq. ([2]) can be 
a weighted average, where each particle is given a weight 
depending on pseudorapidity and/or transverse momen- 
tum (if measured). Our goal here is to characterize 
initial-state fluctuations of the event shape, which are ap- 
proximately independent of rapidity [3] . Weights should 
therefore be chosen independent of (pseudo)rapidity, 
which is a nonstandard choice for odd harmonics [20|. 
With a symmetric detector, one thus selects the rapidity- 
even part of v n . We are concerned with experimental 
observables that can be constructed from v\, i>2 and 113. 
The study of W4 and higher harmonics is more compli- 
cated due to the large interference with V2 2l|, and is 
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left for future work. 

Inserting Eq. ([IJ into Eq. @ gives 

v{ni, . . .,n k } = (v ni . ..v nk cos(ni* ni + . . . + n fe *„J) , 

(3) 

where the average is now only over events. To the extent 
that correlations are induced by collective flow, azimuthal 
correlations measure moments of the flow distribution. 

The simplest v n measurement is the pair correla- 
tion [22}, which corresponds to the event-averaged root- 
mean-square v n 



v n {2} = \AK -n] ~ vW)- 



(4) 



Higher-order correlations yield higher moments of the v n 
distribution: 



v{n,n, -n, -n} = 2v n {2} 4 - v n {4} 4 ~ (w„ 4 ), 



(5) 



where we have used the standard notation u n {4} for the 
4-particle cumulant 23|. 

Finally, one can construct correlations involving mixed 
harmonics, as in previous analyses of i»4 [Hi ] and v\ [24j |. 
The first non-trivial correlations between v\, vi and v 3 
are 



V X 2=v{l,l,-2}, 

v 23 = 42,2,2,-3,-3}, 



v 13 = v{l, 1, 1, -3}, 
v 123 =v{l,2,-3}. (6) 



These observables are new. Note that v{l,l, —2} has 
been analyzed with rapidity-odd weights in harmonic 1[6, 
Il3j |. not with rapidity-even weights. 

One generally expects v\ < v 3 < v%. Thus correlations 
involving high powers of v\ are more difficult to measure. 
As explained in detail in Ref. (2||, the analysis must be 
done in such a way as to isolate the correlation induced 
by collective flow from other "nonflow" effects, which fall 
into two categories: 1) global momentum conservation, 
whose only significant contribution is in v{l, — 1} and 
v{l, 1, —2}. This effect can be suppressed by using the 
p t -dependent weight w = p t — (p+)/(Pt) for at least one of 



the particles in harmonic 1 [18j. 2) Short-range nonflow 
correlations, which can be suppressed by putting rapi dity 
gaps between some of the particles. As shown in Ref. [25], 
all the correlations we have introduced are likely to be 
measurable at the LHC. 



PREDICTIONS 

The anisotropy in the distribution, Eq. (p}, has its ori- 
gin in the spatial anisotropy of the transverse density dis- 
tribution at early times. Following Teaney and Yan fl7j . 
we define 
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FIG. 1. (Color online) Predictions for ratios of various pro- 
posed measurements as a function of centrality (fraction of the 
total cross section, most central to the left) in Au-Au colli- 
sions at RHIC and Pb-Pb collisions at LHC, using a Glauber- 
and a CGC-type model with 100 million and 20 million events, 
respectively. The factors in the denominator are shorthand, 
v n = v n {2}. See text for details. 
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where curly brackets denote an average over the trans- 
verse plane in a single event [26| , weighted by the density 
at midrapidity, and the distribution is centered in each 
event, {re ltp } = 0. In this equation, $„ is the minor ori- 
entation angle (corresponding, e.g., to the minor axis of 
the ellipse for ro=2), and e n the magnitude of the respec- 
tive anisotropy. 

Anisotropic flow scales like the initial anisotropy e n 
and develops along <!>„. It is therefore natural to ex- 
pect that at a given centrality, v n — K n e n and \f„ = 
$„, where K n is a constant that contains all infor- 
mation about the hydrodynamic response to the initial 
anisotropy in harmonic n — in particular, medium prop- 
erties such as the equation of state, viscosity, etc. These 
relations are not exact, but event-by-event hydrodynamic 
calculations have shown that they indeed hold to a good 
approximation for V\ 
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andw 3 [23,1221 — in 
the majority of events, the angles ^ n ~ are strongly 
correlated, while v n /e n is close to a constant value for a 
given set of parameters. There are typically only small 
and apparently random deviations on an event-by-event 
basis, making the use of this approximation very use- 
ful for event-averaged quantities such as those considered 
here. The validity of these relations goes beyond hydro- 
dynamics, and they still hold if the system is far from 
equilibrium [3(J. (On the other hand, these simple rela- 
tions are not valid for higher harmonics such as Vi and 
t>5 [2II, which is why we focus here on n < 3.) 

Inserting these proportionality relations into Eq. (|3|). 
we obtain 

v{m, . . .,n k } = K ni . ..K nk e{ni, . . .,n k }, (8) 

where we have introduced the notation 

e{ni, . . .,n k } = (e ni ...e nh cos(ni$ ni + . . . + n fe $„J) . 

(9) 

Thus the measured correlations are sensitive to details 
of the hydrodynamic evolution mostly through the co- 
efficients K n , and to the initial-state dynamics through 
e{m, . . . , rife}, which only contains information about the 
system prior to equilibration. 

At present, properties of the early-time system are 
poorly constrained by experiment, and contribute the 
largest source of uncertainty in the extraction of medium 
properties like shear viscosity 3lJ. However, these new 
proposed flow measurements will now provide an oppor- 
tunity to significantly constrain the initial state. 

In particular, one can eliminate the dependence on 
the proportionality coefficients K n — and therefore iso- 
late initial-state effects — by measuring the correlations 
defined by Eq. ©, integrated over phase space, and by 



scaling them appropriately: 
u{ni,n 2 , ...,rik} _ 



e{n 1 ,n 2 , . ■ . ,n k } 
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where e n {2} = y/(e^}. The left-hand side of Eq. (JTUJ) can 
be measured experimentally, while the right-hand side 
depends only on early-time dynamics, and can be calcu- 
lated using a model of initial-state fluctuations. Thus, 
although the relations ([5]) are only approximate, taking 
these ratios minimizes sensitivity to medium properties 
and — even though they come from correlations between 
soft particles — they represent some of the most direct 
probes of initial-state dynamics available. 

Eq. (fTOj) holds if the coefficients K n are positive. With 
the definitions Eq. ([7]), this always holds for K 2 and K3. 
However, K\ is negative for low p t particles IH 17 1. One 
can compensate for this negative sign by giving a negative 



18|. 



weight to low p t particles in Eq. @ 

We make predictions for these ratios by computing 
them with two of the most common models for the ini- 
tial state of a heavy-ion collision. First is the PHOBOS 
Glauber Monte-Carlo 



32] 



with binary collision fraction 
x = 0.145 for RHIC collisions and x = 0.18 for LHC. 
The second uses the gluon density from a color-glass- 
condensate (CGC) inspired model — the MC-KLN j^, 
with rcBK unintegrated gluon densities 34|. The main 



difference between the two models is that the eccentric- 
ity is larger in the CGC model 35|, |36j. Both models are 
fairly simple, with the only source of fluctuations being 
the nucleonic structure of nuclei. In reality other sources 
of fluctuations could be important, and future study will 
be needed to fully understand the constraints imposed 
on the initial dynamics by these measurements. 

Fig. [1] displays predictions for all of the scaled corre- 
lations in Au-Au collisions at 200 GeV per nucleon pair 
and Pb-Pb collisions at 2.76 TeV per nucleon pair. 

The top three panels show v{n, n, — n, — n}/v„{2} 4 = 



( v n) / ( v n) 2 ■ For Gaussian fluctuations [3JJ, the n=l and 
3 ratios are equal to 2 (i.e., w„{4}=0), and likewise for 
n—2 in central collisions. However, this is expected only 
in the limit of a large system. A more detailed analy- 
sis [38[ shows that, e.g., t>3{4} should be smaller than 
«3{2} only by a factor ^2 in mid-central collisions, in 
agreement with these results. Note that wherever the 
ratio is greater than 2, the fourth cumulant w n {4} is un- 
defined. The top panel also shows existing data from 
STAR [39] and ALICE 0- Neither measurement in- 
cludes a rapidity gap, and thus may contain nonflow 
correlations (see the discussion in Ref. [HI]). For STAR 
«2{2}, we use both the raw data, and the value with an 
estimated correction for nonflow effects [iol ]. The data 
seem to favor larger relative fluctuations than are con- 
tained in the MC-KLN model used here. 

The bottom four panels display scaled mixed correla- 
tions, indicating non-trivial correlations between ^1, ^2 
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and ^3. The scaled correlation V23 indicates a negligible 
correlation between ^>2 and \&3 up to 40% centrality in 
the Glauber model 0, [4l| , while the CGC model predicts 
a small anticorrelation. In contrast, has both a strong 
correlation with ^3 |42[_f positive Vis) and a (weaker) an- 
ticorrelation with ^2 [121 (negative V12), though this de- 
creases for central collisions. The dependence on impact 
parameter can be attributed to the intrinsic eccentric- 
ity of the nuclear overlap region (jsj . The strong posi- 
tive correlation between ^1 and ^3 explains why v\2 and 
f\<:\ 



17] have the same sign and behave similarly. 



CONCLUSION 



We have proposed a new set of independent flow ob- 
servables in heavy-ion collisions which can be combined 
to tightly constrain theoretical models. In particular, cer- 
tain ratios are constructed which are largely determined 
only by the initial state, and thus directly measure prop- 
erties of the early-time system. We have presented pre- 
dictions for these ratios using two common Monte-Carlo 
models, and compared to existing data. 
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0093-01 and by CEFIPRA under project 4404-2. 



[2] 

[3] 

[4] 

[5] 
[6] 

[7: 

[8 

[9; 

[10 

[11 

[12 



B. Alver and G. Roland, Phys. Rev. C 81, 054905 (2010) 
[Erratum-ibid. C 82, 039903 (2010)] [arXiv:1003.0l94l 
[nucl-th]]. 

B. Alver et al. [PHOBO S Collaboration], Phy s. Rev. 
Lett. 98, 242302 (2007) [arXivmucl ex/0610037] . 

M. Luzum, P hys. Lett. B 696, 499 (2011) 
[arXiv:1011.5773l [nucl-th]] . 

S. Voloshin and Y. Zhang, Z. Phys. C 70, 665 (1996) 
[arXiv:hep-ph/9407282| . 

J.-Y. Ollitrault, Phys. Rev. D 46, 229 (1992). 

C. Alt et al. [NA49 Collaboration], Phys. Rev. C 68, 
034903 (2003) [arXiv:nucl-ex/0303001| . 

K. H. Ackermann et al. [STAR Collaboration], Phys. 

Rev. Lett. 86, 402 (2001) [arXiv:nucl-ex/000901l| . 

S. S. Adler et al. [PHEN IX Collaboration], Ph ys. Rev. 

Lett. 91, 182301 (2003) [arXiv:nucl-ex/0305013| . 

B. B. Back et al. [PHOB OS Collaboration], Phys. Rev. 

C 72, 051901 (2005) [arXiv:nucl-ex/0407012l . 

K. Aamodt et al. [ The ALICE Collaboration ], Phys. 

Rev. Lett. 105, 252302 (2010). |arXiv:1011 .3914 [nucl- 

ex]]. 

A. Adare et al. [PHENIX Collaboration], Phys. Rev. C 
78, 014901 (2008) |arXiv:0801.4545l [nucl-ex]]. 
J. Adams et al. [Star Collaboration], Phys. Rev. C 75, 
034901 (2007) [arXiv:nucl-eI70607003 . 



[13 
[14 

[IS! 

[16 

[17; 

[is; 

[19 

[20 

[21 

[22 
[23 

[24 

[25 

[26 

[27 

[28 

[29 

[30 

[31 

[32 
[33 
[34 
[35 

[36 
[37 
[38 
[39 
[40 

[41 

[42 



J. Adams et al. [STAR Collaboration], Phys. Rev. Lett. 
92, 062301 (2004) [arXiv:nucl-ex/0310029] . 

A. Adare et al. [PHENI X~Collaboration] , Phys. Rev. 
Lett. 105, 062301 (2010) [arXiv: 1003. 55861 [nucl-ex]] . 

F. G. Gardim, F. Grassi, Y. Hama, M. Luzum and 
J.-Y. Ollitrault, Phys. Rev. C 83, 064901 (2011) 
[arXiv: 1103.46051 [nucl-th]]. 

B. B. Back et al. [PHOBOS Colla boration], Ph ys. Rev. 
Lett. 97, 012301 (2006) [arXiv:nucl-ex/0511045| . 

D. Teaney and L. Yan, Phys. Rev. C 83, 064904 (2011) 
[arXiv: 1010. 18761 [nucl-th]]. 

M. Luzum and J.-Y. Ollitrault, Phys. Rev. Lett. 106, 
102301 (2011) arXiv:1011.6361l [nucl-ex]]. 
A. Dumitru, F. Gelis, L. McLerran and R. Venugopalan, 
Nucl. Phys. A 810, 91 (2008) [arXiv: 0804. 38581 [hep-ph]] . 

A. M. Poska nzer and S. A. Voloshin , Phys. Rev. C 58, 
1671 (1998) [arXiv:nucl-ex/9805001] . 

N. Borghini and J.-Y Ollitrau lt, Phys. Lett. B 642, 227 

(2006) [arXiv:nucl-th/0506045| . 

S. Wang et al, Phys. Rev. C 44, 1091 (1991). 

N. Borghini, P. M. Di nh and J.-Y. O llitrault, Phys. Rev. 

C 64, 054901 (2001) [arXiv:nucl-th/0105040| . 

N. Borghini, P. M. D inh and J.-Y. Ollitrault ^ Phys. Rev. 

C 66, 014905 (2002) [arXiv:nucl-th/0204017| . 

R. S. Bhalerao, M. Luzum and J. Y. Ollitrault, J. Phys. 

G 38, 124055 (2011) |arXiv:1106.4940l [nucl-ex]]. 

R. S. Bhal erao and J.-Y. Ollitrault , Phys. Lett. B 641, 

260 (2006) [arXiv:nucl-th/0607009j . 

H. Holopainen, H. Niemi and K. J. Eskola, Phys. Rev. C 
83, 034901 (2011) arXiv: 1007.03681 [hep-ph]]. 
Z. Qiu and U. W. Heinz, Phys. Rev. C 84, 024911 (2011) 
[arXiv:1104.0650l [nucl-th]]. 

G. Y. Qin, H. Petersen, S. A. Bass and B. M uller, Phys. 
Rev. C 82, 064903 (2010) [arXiv:1009.1847l [nucl-th]].' 
N. Borghini and C. Gom beaud, Eur. Phys. J. C 71, 1612 
(2011) |arXiv:1012.0899l [nucl-th]]. 

M. Luzum and P. Romatschke, Phys. Rev. C 78, 
034915 (2008) [ Erratum-ibid. C 79, 039903 (2009)] 
[arXiv:0804.40T5l [nucl-th]]. 

B. Alver, M. Baker, C. Loizides and P. Steinberg, 
larXiv:0805.44"TD [nucl-ex]. 

H. J. Drescher and Y. N ara, Phys. Rev. C 76, 041903 

(2007) [arXiv:0707.0249l [nucl-th]]. 

J. L. Albacete and A. Dumitru, arXiv:1011.5161 [hep- 
ph]. 

T. Hirano, U. W. Heinz, D. Kharzeev, R. Lacey 
and Y. Nara, Phys . Lett. B 636, 299 (2006) 
[arXiv:nucl-th /051 1046] . 

T. Lap pi and R. Ven ugopalan, Phys. Rev. C 74, 054905 
(2006) [arXiv:nucl-th/0609021| . 

S. A. Voloshin, A. M. Poskanzer, A. Tang and G. Wang, 
Phys. Lett. B 659, 537 (2008) [arXiv:070 8.0800 [nucl-th]]. 
R. S. Bhalerao, M. Luzum and J. -Y. Ollitrault, Phys. 
Rev. C 84, 054901 (2011) [arXiv:1107.5485l [nucl-th]]. 
J. Adams et al. [STAR C ollaboration] Phys. Rev. C 72, 
014904 (2005) [arXiv:nucl-ex/0409033| . 
J.-Y. Ollitrault, A. M. Poskanzer and S. A. Voloshin, 
Phys. Rev. C 80, 014904 (2009) |arXiv:0904.23T5l [nucl- 
ex]]. 

J. L. Nagle, M. P. McCu mber, Phys. Rev. C83, 044908 
(2011). |arXiv:1011. 18531 [nucl-ex]]. 

P. Staig and E. Sh uryak, Phys. Rev. C 84, 034908 (2011) 
[arXivil 008.31391 [nucl-th]]. 



